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Abstract 

The purpose of this paper is to determine quantum master and filter equations for 



C^ ' systems coupled to fields in certain non-classical continuous-mode states. Specifically, 



we consider two types of field states (i) single photon states, and (ii) superpositions 
of coherent states. The system and field are described using a quantum stochastic 
unitary model. Master equations are derived from this model and are given in terms 
^ ■ of systems of coupled equations. The output field carries information about the 

'nI" , system, and is continuously monitored. The quantum filters are determined with the 

aid of an embedding of the system into a larger non-Markovian system, and are given 

ps| ■ by a system of coupled stochastic differential equations. 

I> 

O I Keywords: quantum filtering, continuous-mode single photon states, continuous-mode 

superpositions of coherent states, quantum stochastic processes 

^ '. 1 Introduction 



In recent years single photon states of light and superpositions of coherent states have 
become increasingly important due to applications in quantum technology, in particular, 
quantum computing and quantum information systems, [25], [27], [23], [18], [30]. For 
instance, the light may interact with a system, say an atom, quantum dot, or cavity, and 
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Single Photon Quantum Filter 2 

this system may be used as a quantum memory, [25], or to control the pulse shape of 
the single photon state [27]. When light interacts with a quantum system, information 
about the system is contained in the scattered light. This information may be useful for 
monitoring the behavior of the system, or for controlling it. The topic of this paper concerns 
the extraction of information from the scattered light when the incoming light is placed in 
a single photon state |\I^) = |1^), or a superposition of coherent states |\1/) = ^ ■ aj\fj), as 
illustrated in Figure [H 
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Figure 1: A system initialized in a state \rf) coupled to a field in a state |\E') (single photon or 
superposition of coherent states). The output field is continuously monitored by homodyne 
detection (assumed perfect) to produce a classical measurement signal Y{t). The output 
Y{t) is filtered to produce estimates X{t) = 7r((X) of system operators X at time t. 

The problem of extracting information from continuous measurement of the scattered 
light is a problem of quantum filtering, [1], [5], [6], [H], [31], [3], [9], [32]. The current state 
of the art for quantum filtering considers incoming light in a vacuum or other Gaussian 
state, with quadrature or counting measurements. Both single photon states of light, and 
superpositions of coherent states of light, are highly non-classical, and are fundamentally 
different from Gaussian states. In view of the increasing importance of these non-Gaussian 
states of light, the purpose of this paper is to solve a quantum filtering problem for systems 
driven by fields in single photon states and superpositions of coherent states. 

In the case of single photon fields, the master equation describing unconditional dy- 
namics was shown to be a system of coupled equations in [17], a feature of non-Markovian 
character. Markovian embeddings were used in [TO] to derive quantum trajectory equa- 
tions (quantum filtering equations) for a class of non-Markovian master equations. In recent 
work, the authors have shown how to construct ancilla systems to combine with the system 
of interest to form a Markovian extended system driven by vacuum from which quantum 
filtering results may be obtained for single photon states and superpositions of coherent 
states from the standard filter for the extended system, p^, [20]. However, depending on 
the complexity of the non-classical state, it may possibly be difficult to determine suitable 
ancilla systems, and indeed the superposition case was not straightforward. In this paper 
we present an alternative approach to the embedding that also allows for the derivation of 
the quantum filter. The extended system forms a non-Markovian system, with the ancilla, 
system and field initialized in a superposition state. While standard filtering results do 
not apply, the quantum stochastic methods can nevertheless be applied to determine the 
quantum filters. In this way, we expand the range of methods that may be applied to 
derive quantum filters for non-classical states. 
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The paper is organized as follows. In section [2] the idealized filtering problem to be 
solved in this paper is formulated. The continuous mode single photon states are defined 
in Section 13.11 The master equation for the single photon field state is derived in Section 
13.21 using the model presented in Section El This leads naturally to Section 13.31 where the 
system is embedded in a larger model, inspired by the approach used in [TO] but differing in 
the details. This extended system provides a compact and transparent description of the 
problem, and may readily be generalized to n-photon states, and indeed, multiple channels 
of n-photon states. The quantum filter for the extended system is presented in Section [3^ 
with a derivation extending the reference method appearing in Section 13.61 The filtering 
results for the extended system are used to find the filtering equations for the original 
problem involving a single photon field in Section 13.51 The superposition of coherent field 
states is defined in Section 14. H and a suitable embedded system for this case is described 
in Section 14.21 The corresponding master and filtering equations are presented in Sections 
14.31 and 14. 4[ respectively. Some concluding remarks are made in Section O 

In this paper we are not concerned with technical issues concerning domains of un- 
bounded operators and related matters, and indeed, we assume that the system operators 
are bounded, and that all quantum stochastic integrals are well-defined in the sense of 
Hudson-Parthasarathy, |22j . 

Notation: We use the standard Dirac notation \tp) to denote state vectors (vectors in 
a Hilbert space) [26], [1]. The superscript * indicates Hilbert space adjoint or complex 
conjugate. The inner product of state vectors \ipi) and \ip2) is denoted {4'i\ip2)- The 
expected value of an operator X when the system is in state 1-0) is denoted ]E^[X] = 
{i/j\X\i/)). For operators A and B we write {A,B) = tT[A*B]. 

2 Problem Formulation 

We consider a quantum system 5* coupled to a quantum field B, as shown in Figure [H The 
field B has two components, the input field Bin and the output field. Bout- In this paper we 
consider two non-classical cases for the state |^) of the input field (i) a single photon state 
1^) = |1^), where (^ is a complex valued function such that J^ \^{s)\'^ds = 1 (representing 
the wave packet shape), or (ii) a superposition of coherent states I^P) = J2j'^j\fj)y where 
\fj) are coherent states and the complex numbers aj {j = 1, . . . ,n) are normalized weights. 

As illustrated in Figure [H the field interacts with the quantum system S, and the 
results of this interaction provide information about the system that may be obtained 
through continuous measurement of an observable Y{t) of the output field Boutif)- The 
filtering problem of interest in this paper is to determine the conditional state from which 
estimates X{t) of system operators X may be determined at time t based on knowledge of 
the observables {Y{s), < s < t}. 

In what follows the system S is assumed to be defined on a Hilbert space S^s^ with an 
initial state denoted \rf) G S)s- The input field Sj„ is described in terms of annihilation 
B[^) and creation B*{^) operators defined on a Fock space 5^, [2Sl Chapter II], [HI Section 
4]. Quantum expectation will be denoted by the symbol E, and when we wish to display 
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the underlying state, we employ subscripts; for example, E^,j, denotes quantum expectation 
with respect to the state \r]) ® |\l/). 

The dynamics of the system will be described using the quantum stochastic calculus, 
[22] . [T5] . [28] . [16], [9]. Quantum stochastic integrals are defined in terms of fundamental 
field operators B{t), B*{t) and A(t), f231 Chapter II], ^ Section 4]|ll The non-zero Ito 
products for the field operators are 

dB(t)dB*{t) = dt, dB{t)dA{t) = dB{t), dA{t)dA{t) = dA{t), dA{t)dB*{t) = dB*{t). 

(1) 
The dynamics of the composite system is described by a unitary U{t) solving the 
Schrodinger equation, or quantum stochastic differential equation (QSDE), 

dU{t) = {{S - I)dA{t) + LdB*{t) - L*SdB{t) - {-L*L + iH)dt}U{t), (2) 

with initial condition ^(0) = /. Here, H is a fixed self-adjoint operator representing the 
free Hamiltonian of the system, and L and 5* are system operators determining the coupling 
of the system to the field, with S unitary. In this paper, for simplicity we assume that 
the parameters S, L, H are bounded operators on the system Hilbert space S^s- However, 
we remark under some suitable additional conditions the results and equations obtained 
in this paper should also be extendable to some special classes of QSDEs with unbounded 
parameters, exploiting the results in [131 HI] • 

A system operator X at time t is given in the Heisenberg picture by X{t) = it{X) = 
U{t)*{X I)U{t) and it follows from the quantum Ito calculus that 

d3t{X) = US*XS-X)dAit)+jtiS*[X,L])dBity 

+jti[L*,X]S)dB{t)+jtiCiX))dt, (3) 

where 

CiX) = h*[X, L] + i[L*, X]L - ^[X, H]. (4) 

The map X i— )■ C{X) is known as the Lindhlad generator, while the quartet of maps 
X H- £(X), S*XS — X, S*[X, L], [L*, X]S are known as Evans-Hudson maps. 

The output field is defined by Bout{t) = U{t)*B{t)U{t)u In this paper we consider the 
output field observable Y{t) defined by 

Y{t) = U{tYZ{t)U{t), (5) 

where 

Z{t) = B{t) + B*{t), (6) 



""^In terms of annihilation and creation white noise operators 6(i), h*{t) that satisfy singular commutation 
relations [h{s),h*{t)\ — 5{t — s), the fundamental field operators are given by B{t) = J^ b{s)ds, B*{t) = 
/q h*{s)ds, and A(t) = /J b*{s)b{s)ds. Also, we may write B{() = J^ Cis)dB{s). 

^Recall B{t) — -Bm(t) is the input field. 
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is a quadrature observable of the input field (the counting case Z(t) = A(t) is discussed 
briefly in Section [5]). Note that both Z{t) and Y{t) are self-adjoint and self-commutative: 
[Z{t), Z{s)] = and [Y(t), Y{s)] = 0. We write ^ and ^ for the subspaces of commuting 
operators generated by the observables Z{s), Y{s), < s < t, respectivelyjfl They are 
related by the unitary rotation ^ = U [t]* S'tU {t) . Physically, Y{t) may represent the 
integrated photocurrent arising in an idealized (perfect) homodyne photodetection scheme, 
as in Figure [H For further information on homodyne detection, we refer the reader to the 
literature; for example, [2], [3], [32] . 

The primary goal of this paper is to determine the quantum filter for the quantum 
conditional expectation (see, e.g. P, Definition 3.13]) 

l(t)=E,*[X(t)|^]. (7) 

This conditional expectation is well defined, since X{t) commutes with the subspace ^ 
(non-demolition condition). The conditional estimate X{t) is affiliated to ^ (written in 
abbreviated fashion as X{t) G ^) and is characterized by the requirement that 

E^vi, [X{t)K] = Er,q, [X{t)K] (8) 

for all K e^f. 



3 Single Photon Input Fields 

3.1 Single Photon Fields States 

In this section we consider the continuous-mode single photon state |\1/) = |1^) defined by 
[211 sec. 6.3], [271 eq. (9)] 

\U)=B*mO), (9) 

where (^ is a complex valued function such that f^ \^{s)\'^ds = 1, and |0) is the vacuum 
state of the field. Expression ([9]) says that the single photon wavepacket with temporal 
shape ^ is created from the vacuum using the field operator B*{C,). 

The Hilbert space for the composite system is 

where here we have exhibited the continuous temporal tensor product decomposition of 
the Fock space d = dt]^ d{t into past and future components, which is of basic importance 
in what follows. Write 

Eu[X®F] = {r^miX ^ F)\7^^) = (r7|X|r/)(l^|F|l^) (10) 

for the expectation with respect to the product state jr/l^), where the field is in the single 
photon state. Here and in what follows X is a bounded system operator acting on S^s, and 



"^^ft and ^ are commutative von Neumann algebras. They are also filtrations, e.g. ^^ C ^2 whenever 
ti < t2. 
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F is a field operator acting on the Fock space 5- Similarly, we may define the expectation 
when the field is in the vacuum state, 

Eoo[X® F] = {7^0\{X ® F)\7]0) = {r]\X\r]){0\F\0). (11) 

We will also have need for the cross-expectations 

Eio[X ®F] = (r/l^|(X ® F)\r]0), and Eoi[X ® F] = (r/0|(X ® F)\7]l^). (12) 

A crucial difference between the single photon state and the vacuum state is that the 
later state factorizes |0) = |0i])(8)|0(i) with respect to the temporal factorization 5" = dt]'^d{t 
of the Fock space, with |0i]) G 5t] and |0(() G ^(t, while the former does not. Rather, we 
have 

\k) = B*{O\0) = 1^) ® |0(,) + |0,]) ® |1^(J, (13) 

where 

\ht]) = BrmOt]), and |l5(,) = Bt*mOit), (14) 

and 

BiiO = BiCxio,t])^ BtiO = Bi^X(t,oo]), BiO = BiiO + BtiO- (15) 

Here, X[Q,t] is the indicator function for the time interval [0,t]. Note that while |1^) has 
unit norm, we have 

II |lc*]) f = / \as)?ds, and II |1^(,) f = / \as)\'ds. (16) 

Jo Jt 

A consequence of the additive decomposition ( IT3l) and the definitions ( TT5|) is the follow- 
ing. Let -ft'(t) be a bounded operator acting on the full Hilbert space i^ that is adapted, 
i.e. K{t) acts trivially on ^{ti the field in the future. Then the expectation with respect to 
the single photon field may be expressed in terms of the vacuum state as follows: 

En[/r(t)] = ¥.oo[Bi{i)K{t)Bm)+r{t)K{t)] (17) 

where r{t) = f^ \^{s)\'^ds. 

3.2 Master Equation 

Before deriving the quantum filter, we work out dynamical equations for the unconditioned 
single photon expectation, [17]. To assist us in evaluating this expectation, we make use 
of the following lemma. 

Lemma 3.1 Let K{t) be a bounded quantum stochastic process defined by 

K{t) = f Mo{s)ds+ f M_{s)dB{s)+ f M+{s)dB*{s)+ f M^{s)dA{s), (18) 
Jo Jo Jo Jo 
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where Mq, M± and Mi are bounded and adapted. Then we have 

Kn[K{t)] = En[/ Mo{s)ds] +Eio[ [ M_(s)e(s)ds] 
Jo Jo 

+Eoi[/ M+{s)C{s)ds]+Eoo[f Mi{sms)\^ds], (19) 

Jo Jo 

Eio[K(t)] = Eio[ / Mo{s)ds] +Eoo[ I M+{s)C{s)dsl (20) 

Jo Jo 

Eoi[K(t)] = Eoi[ / Mo(s)rfs] +Eoo[ / M.{s)i{s)dsl (21) 

Jo Jo 

Eoo[K{t)] = Eool [ Mo{s)ds]. (22) 

Proof. Using ([IT]), the expressions 5^(0 = IoCis)dB{s), B^*{^) = J^Cis)dB*{s), 
and the Ito rule we have 

Eu[dK{t)] = Eoo[diBi{OKit)B-*iO + r{t)Km 
= Eoo[Bi{OMo{t)B-*{0 + r{t)Mo{t) 

+M+it)B-*m*it) + BiiOM-it)m + Miitmt)\']dt 

= En[Mo{t)]dt + Eoo[M+{t)B*{OC{t) + B{OM_{t)m + Mi\m\']dti23) 
This last line is justified since M± are adapted and Eoo[-Bt^(0] = 0- That is, 

En[dK{t)] = En[Mo{t)]dt + Eoi[M+{t)]C{t)dt + Eio[M4tMt)dt + Eoo[ [ Mi{sms)\^ds]. 

Jo 

(24) 

This proves ( !T9|) . The remaining expressions are proven in a similar manner. D 

We will first express the master equation in Heisenberg form using the expectations 

^4\X)=E,,[X{t)]. (25) 

Note that for all t > we have 

f^TiI) = l = f^l\l), f^ni) = = f^lV)- (26) 

Theorem 3.2 The master equation in Heisenberg form for the system when the field is in 
the single photon state |lg) is given by the system of equations 

+f,T{S*XS~Smt)\', (27) 

Af (X) = ^il\C{X)) + fi^iS^lX, L])C{t), (28) 

A°^(X) = ^°i(£(X)) + /ir([^%^]5)e(t), (29) 

fiTix) = ^r(^(^))- (30) 

The initial conditions are 

/i^i(X) = ^°°(X) = (r7,Xr/), ^J°(X)=^°i(X)=0. (31) 
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Equations ([2 



( 130|) are obtained by applying Lemma [3. II to the Heisenberg 

D 



Proof. 
equation (^ 

It is apparent from Theorem 13.21 that the single photon expectation fij^{X) = Eii[X(t)] 
cannot be determined by a single differential equation, and that instead a system of coupled 
equations is required, equation (P7|) - (I5U]) . Note that the unitary matrix S appearing in the 
Schrodinger equation ([2]) does appear in the single photon master equations fl7r|) - (15nil . in 
contrast to the vacuum case (which corresponds to fl30l) ). 



3.3 Embedding 

In this section we construct a suitable embedding for the system and single photon field, and 
show how the system of master equations from Section [3l2] can be compactly represented as 
a single equation for a larger system. This embedding will be used in subsequent sections 
to derive the quantum filter. We should emphasize, however, that our embedding is not 
the same as that used in [10], [19], [20]. The embedding is illustrated in Figure [2l 
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Figure 2: System embedded in the extended system. While the analysis does not employ 
any coupling between the system and ancilla two-level system, the ancilla, system and field 
are assumed to be initialized in a superposition state |S) defined in equation fl36|) . 

Recall that the system and field are defined on a Hilbert space i^ = i35 (g) 5- We define 
an extended space 

S^ = C'(^S)=S)®S) (32) 

which includes the system, field and an ancilla two-level system. Let |eo) and |ei) be an 
orthonormal basis for C^, 



|eo) 



|ei) 



and let A be an operator acting on C^, i.e. a complex 2x2 matrix 

. _ Oil ttio 



(33) 



(34) 
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It may be helpful to think of operators A0 X F on the extended space S) represented in 
the Kronecker product form 



A®{X®F) 



aiiiX^F) aioiX^F) 
aoi{X®F) aoo{X ®F) 



(35) 



We allow the extended system to evolve unitarily according to J ® U{t), where U{t) 
is the unitary operator for the system and field, given by the Schrodinger equation ([2]). 
Note in particular that the system is not coupled to the ancilla C^, and observables of this 
two-level system are static. We initialize the extended system in the superposition state 



|S) = aildr^l^) + ao|eo^O), 
where |aoP + |«ip = 1. This state evolves according to 

|S(t)) = (/®t/(t))|S). 
For notational convenience we write 



and note that w = J2jk'^jk\^j)(^k\ is a density matrix for C^. 

The expectation with respect to the superposition state |S) is given by 

j2t{A®X) = E^[A®X{t)] = (S|(A®X(t))|E) = ^ W7,,a,fe^f (X). 

jk 



(36) 



(37) 



(38) 



(39) 



This expectation is correctly normalized, /i((/ (X> /) = 1, and the expectations jjI (X) 
defined in Section [3^2] are scaled components of /if (A X): 



Jk 



f^nx) 



fj^t{\ej){ek\ ^X) 



Wjk 



(40) 



for Wjk 7^ 0, otherwise it can be set to, say, 0. We also have 



jfc 



/^r(^) 



wufiti\ej){ek\ OX) 
Wjkfiti\ei){ei\ ®/) ' 



(41) 



Note that in the extended space the Schrodinger and Heisenberg pictures are related 



by 



EE(t) [A®X®F]=Ej:[A®U* (t) (X ® F) [/(t)] . 



(42) 



In order to derive the equation for expectations in the extended system, we need the 



following lemma, which follows from Lemma I3.H and makes use of the matrices 



f^+ 



|ei)(eo| 



1 




a^ 



|eo)(ei| 





1 



(43) 
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Lemma 3.3 Assume ao 7^ 0, and let M(t) he hounded and adapted. Then 

Es[ / A ® M{s)dB{s)] = z/Es[ / (Aa+) ® M{s)^{s)dsl (44) 

Jo Jo 

E^lf A^ M{s)dB*{s)] = z/*Es[ / (a_A) ® M{s)C{s)ds], (45) 

Jo Jo 

Es[ / A ® M{s)dA{s)] = |z/|2Ee[ /" (a_A(7+) ® M(s)|e(s)prfs], (46) 



where 



Qio' 



(47) 



In Lemma 13.31 expectations of stochastic integrals with respect to the superposition 
state |E) are expressed in terms of expectations of non-stochastic integrals again with 
respect to |S) with the aid of the matrices a± acting on the ancilla system C^. The action 
of the field annihilation, creation and gauge operators is therefore captured algebraically 
and all expectations in these relations are with respect to the same state. 

We now have 

Theorem 3.4 Assume ao 7^ 0. Then the expectation jlt{A® X) (defined hy (3^) ) evolves 
according to 

iit{A(g)X) = jlt{gt{A^X)), (48) 

where 

GtiA^X) = A®C{X) + {Aa+)®[L*,X]Su^{t) + {a^A)®S*[X,L]u*C{t) 

+ ia_Aa+)0iS*XS-X))\u^it)\\ (49) 

The reader may easily verify that the system of master equations (l27 1) -( l30l) for /i^ {X), 
j, k = 1,0, follows from equation fH5]) by setting A = \ej){ek\. 

3.4 Quantum Filter for the Extended System 

The extended system provides a convenient framework for quantum filtering, since all 
expectations can be expressed in terms of the superposition state |S). Our immediate goal 
in this section is to determine the equation for the quantum conditional expectation 

nt{A^X)=Ej,[A^X{t)\I^^t], (50) 

and in Section 13.51 we will explain how the quantum filter for the single photon field may 
be obtained from this equation. 

The continuously monitored field observable that corresponds to the conditional expec- 
tation ( l50|) is / ® y{t), and from ^ we have the corresponding output equation for the 
extended system: 

d{I Y{t)) = / ® (L(t) + L*{t))dt + 1® {S{t)dB{t) + S*{t)dB*{t)). (51) 
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In what follows we will make use of the following lemma concerning expectations of the 
process 

V{t) = f {S{s)dB{s) + S*is)dB*{s)) (52) 

Jo 

with respect to the single photon state. 
Lemma 3.5 For any K E 'S^s, we have 

En[{V{t) - V{s))K] = Eio[ / S{r)^{r)dr K] +Eoi[ f S*{r)C{r)drK]. (53) 

J s J s 



Proof. Equation (!53l) is obtained using the additive decomposition ( !T3ll and Lemma 

D 

Note that an operator K in the unital commutative algebra I ^ '3^t has the form K = 

1 ^ K, where K e ^. By the spectral theorem, [9l Theorem 3.3], we may identify K 
and K, both of which are equivalent to a classical stochastic process Kt{s), < s < t. In 
the remainder of this paper, we use these identifications without further comment. The 
quantum conditional expectation TTt{A(^X) G J®^ is well defined because A^X{t) is in 
the commutant / ® ^' of the algebra I ®'3^f, and is characterized by the requirement that 

Es[^t(A ® X)I ® K] = Es[(A X{t)){I ® K)] (54) 

for all K e^t, see, e.g. p Definition 3.13]. 

Theorem 3.6 Assume ao ^ 0. The conditional expectation 7Tt{A^X) defined by (E^ for 
the extended system satisfies 

dnt{A ®X) = ntiGtiA ® X))dt + UtiA ® X)dW{t), (55) 

where 

UtiA ® X) = ^t{A ® {XL + L*X)) - nt{A ® X)TVt{I ® (L + L*)) 
+^i((A(7+) ® XS)u^{t) + 7fi((a_A) ® S*X)iy*C{t) 
-^t{A ® X)^i((a+ 8) 5)z/e(t) + (a_ ® S*)u*C{t)) (56) 

anJ 

diy(t) = dY{t) - nt{I O (L + L*) + (a+ ® 5)i^^(t) + (a_ ® 5*)//*r(t))dt. (57) 



r/ie process W{t) defined by ( [57| j is a I ® 'S^t Wiener process with respect to |S) anc? zs 
ca/Zed the innovations process. 

Proof. We follow the characteristic function method [29], [7], [5], whereby we pos- 
tulate that the filter has the form 

d^t{A ®X)= 7t{A ® X)dt + Ut^A ® X)I ® dY{t), (58) 
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where J-'t and Tit are to be determined. 

Let / be square integrable, and define a process c/ by 

dcfit) = f{t)cf{t)dY{t), Cf{0) = 1. (59) 

Then / (g) c/(t) is adapted to / ® ^, and the defining relation (IMl) imphes that 

Ej,[A ® (X(t)c/(t))] = EE[ft(A ® X)/ ® c/(t))] (60) 

holds for all /. By calculating the differentials of both sides, taking expectations and 
conditioning we obtain 

Ej,[A ® {dX{t)cfit))] = Es[(/ ® Cfit))^t{giA (S) X)) (61) 

+ (/ ® fit)cf{t)){^t{A ® (XL + L*X)) 
+^i(Aa+ ® XS)uC{t) + ^i(a_A (g) S*X)u*C{t)}]dt 

and 

Es[A®(d^i(A®X)c/(t))] (62) 

= Es[(/ ® c/(t){J-i(A ® X) + ?{t(A ® X)^t(/ ® (L + L*)) 
+Hi(A ® X)7ri((a+ ® 5)i/e(t) + (a_ ® 5*)z/*r(i))} 
+ (/ ® fit)cf{t)){ntiA ® X)^i(/ ® (L + L*)) + -Ht 
+nt{A X)^i(a+ ® 5z/e(0 + a_ ® 5V*r(i))}]^^- 

Now equating coefficients of c/(t) and f{t)cf{t) we solve for J-i(A ® X) and ^(^ ® X) to 
obtain the filter equation. 

We now prove the martingale property Es[/ ® {W{t) — W{s)) \ I ® '3^s] = 0, that is, 
E2[/ ® {W{t) - W{s)){I ®K)] = Q for all K e ^f Now 

E^[I®{W{t)-W{s)){I®K)] 
= E^[{I®{Y{t)-Y{s)) 

- I Tiril ® (L + L*) + a+® S^{r) + a_ ® S*Cir))dr}I ® K] 

= Ej:[{I®{Y{t)-Y{s)) 
t 
(/ ® (L(r) + L*{r)) + a+ ® 5(r)e(r) + (T_ ® ^*(r)r(r))rfr}/ ® K] 



= Es[{/ ® {V{t) - V{s)) - I {a+® S{r)p^{r) + a_ ® S*{r)p''C{r))dr]I ®K] = Q. 

To see that this last expression is zero, we make use of Lemma 13. 5^ the multiplicative 
factorization of the vacuum state, the fact that V{t) has zero expectation in the vacuum 
state, to find that 

Es[/ ® {V{t) - V{s))I ®K]= wnEn[(y(t) - V{s))K] + w;ooEoo[(^(t) - V{s))K] 

= w;n( I Eio[5(r)K]e(r)dr+ / E^^[S*{r)K]C{r)dr), 
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and 

Es[ / (a+ ® Sir)K)^{t) + a_ ® S*ir)KCit))dr] 

J s 

= wnif E,o[Sir)K]ar)dr+ f Eqi [S*(r)K]r(r)dr). 



Finally, since dW{t)dW{t) = dt, Levy's Theorem implies that W{t) is a ^ Wiener process. 
This completes the proof. D 

Notice the terms involving a± in the filter (equation flS^ ) and in the innovations process 
(equation flFTI) ). These terms arise from expectations involving the single photon state. 
Note that due to the martingale property of the innovations process W{t) we see that 
if we take the expected value of equation ( 156|) we recover equation ( HHj) . consistent with 
Es[^j(yl ® X)] = flt{A^ X) and the definition of conditional expectation. 

3.5 Single Photon Quantum Filter 

We return now to the main goal of the paper, namely the determination of the quantum 
filter for the conditional state when the field is in the single photon state, as stated in 
equation ([7]). As discussed earlier, our strategy is to make use of the filtering results 
obtained in Section 13.41 for the extended system. 

Lemma 3.7 Assume a^ ^ 0. Define the conditional quantities nl {X) by 

ikfy. wii^t{\ej){ek\®X) 
WjkT^t[\ei){ei\ ®I) 

where Trt{A®X) is the conditional state for the extended system defined by [5^) . Then for 
all K E '3^t we have 

Eu[ntiX)K] = E,k[jtiX)K]. (64) 

Proof. We have 

En[nl\X)K] = —E^[\e,) {e^\ ^ {7rl\X) K)] 

Wii 

= —E^[7rti\e^){ei\ ® /)(/ ® vrf (X)K)] 

Wii 

= — EsK(|e,)(e,|®X)(/®K)] 

Wjk 

= —Ej,[{\e,){ek\®jt{X)){I®K)] 

= E,u[3t{X)K] 

as required. D 

We can now present our main theorem for the quantum filter for the single photon field 
state. 
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Theorem 3.8 The quantum filter for the conditional expectation with respect to the single 
photon field is given in the Heisenberg picture by 

X{t)=Eu[Xit)\Wt]=nl\X), (65) 

where n}^ {X) is defined by (E^) (for j = k = 1), and is given by the system of equations 

dnl\X) = inl\CiX)) + n^\S*[X,L])Cit) + nl'{[L\X]S)m 

+7r^\S*XS-X)\^{t)\^)dt 

H-^l'iXL + L*X) + 'Kf{S*X)C{t) + Til\XS)i{t) 

-nl\X)inl\L + L*) + n^\S)m + 7rl'iS*)CmdWit), (66) 

dTrl^X) = {nl'{C{X)) + n^\S*[X,L])C{t))dt 

+ (7rf (XL + L*X) + 7c'^%S*X)C{t) 

-7rl\X){nl\L + L*) + 7r^\S)m + vrf (5*)r(t)))dH^(t), (67) 

dn^^X) = {n^\C{X)) + n',\[L\X]S)m)dt 

+ {n^\XL + L*X) + 7T^''{XS)m 

-4\X){nl\L + L*) + n^\S)m + nl\S*)CmdW{t), (68) 

dn^^X) = Txf'{C{X))dt+{4%XL + L*X) 

-n^\X){nl\L + L*) + 7r^\S)at) + nl\S*)C{t)))dW{t). (69) 

Here, the innovations process W{t) is a'3/t Wiener process with respect to the single photon 
state and is defined by 

dWit) = dY{t) - {7rl\L + L*) + nl'{S)m + n^\S*)C{t))dt. (70) 

The initial conditions are 

nl\X) = 4'{X) = {v, Xv), 4'{X) = 7r°^(X) = 0. (71) 

Proof. Suppose first that ao ^ 0. Setting j = k = 1 in equation flMl) above, 
and noting that K E '3^s was otherwise arbitrary, we deduce that nl^^X) is the desired 
conditional expectation for the single photon field state, as characterized by equation ([8]). 
The differential equations (!66l) - (!69l) follow from the definition (163|) . the filter (!55l) for the 
extended system, and the Ito rule. Next, we note that the coefficients of the QSDEs (!66l) - 
(169|) . the initial conditions, and Yf do not depend on oq and ai. Hence, the solutions 
Ti-j. (X) of this system of equations are independent of ao and ai. Therefore, vr^ (X) can 
be defined for aj G {0, 1}, j = 0, 1, and is in fact identical for all < |ao|) |oi| < 1- 
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We now prove that W{t) is a ^-martingale, that is, Eii[W{t) - W{s) |^] = 0. To this 
end, let K e ^s- Then 

^ii[{W{t) - W{s)}K] = Eu[{Y{t) - Y{s) 

^\nl\iL + L*) + ^l\S)m + iTf{S*)C{t))dr}K] 



En[{ j {L{r) + L*{r))dr + ¥(1) - V(s) 

J s 

- fi-^l'HL + L*) + nl'{l)m + <\l)C{t))dr}K] 

J s 

EnKy {L{r) + L*{r) - ti]^{{L + L*))dr}K] 

+En[{F(t) - V{s) - f{nl\S)m + n'^\S*)C{t))dr}K] 

J s 



however, this vanishes from (IMI) . Lemma [3 .St and 



= Eio[ f S{r)i{r)drK] + Eoi[ /" 5*(r)r(r)drK]. 



Finally, since dW{t)dW{t) = dt, Levy's Theorem implies that W{t) is a ^ Wiener 
process. D 

3.6 Reference Method for Filtering in the Extended System 

The reference method is one of the standard approaches to filtering theory, with its origins 
in the work of Duncan, Mortensen, Zakai, Holevo and Belavkin, see [12], [21], [5], [8], [9]. In 
this section we apply this approach, as described in [9l sec. 6], to the filtering problem in the 
extended system, giving an independent derivation of the fundamental filtering equation. 

Our first step is the following. 

Lemma 3.9 Assume ckq 7^ 0. Then we have 

Es[(A ® X{t))] = Ej,[F*{t){A ® X)F{t)] (72) 

where F{t) G (/ ^t)' is given by 

dF{t) = {Go{t)dt + Gi{t)dZ{t))F{t), (73) 

F(0) = /, and 

Go{t) = -I® {^L*L + iH) -a+®{L + L*S)v^{t), (74) 

Giit) = I®L + a+®{S- I)ui{t). (75) 
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Proof. Let us suppose that F (t) satisfies ( 173]) with coefficients Go (i), Gi (t) which 
both commute with (T+ ® /, that is, 

G, (t) = / Qio (t) + a+ ® gn (t) . 

Then 

Es [d {F (t)* A ® XF (t)}] = Es [F (t)* {Tt(A ® X)dt + gi(A ® X)rfZ(t)} F (t)] , 

where 

Ti (A ® X) = ^(t)*^ ® XG'i(t) + A XG'o(t) + G'o(t)M ® X, (76) 

Qi (A ® X) = A ® XG'i(t) + Gi{t)*A ® X. (77) 

Using lemma 13.31 we see that this equals 

Ej,[F {ty {Tt{A ® X) + Qt{A ® X)u^ (t) (a+ ® /) + i/*^ (0* (^- ® /) Qt(^ ® X)}F {t)]dt. 

We now require that this equals Es [^t {A (g) X)] (it. This implies the four identities 

gn*Xg^, + y^gl^X + y^CXg^i = HWS*XS - X) (78) 

gl^Xg^^ + X(7oi + ^^^^1*0^ + i^^Xg^^ = v^[L\X]S (79) 

^*iX^io + goi*X + z/*r^io*X + v*CXgio = v*CS*\X, L] (80) 

gl^Xg,, + ^*o^ + ^^00 = C{X). (81) 

The first identity flTHl) is satisfied if (711 = /^(^(S* — 1). Substituting X = I into fl79|) we 
deduce that goi = —^iglQS — v^gw and thus 

9wXgii + XgQi + z/^s-ioX + z/^X^io = z^^bio, ^l-S*. 

Therefore (!79|) is satisfied if g^ = L, and consequently goi = —i/^{L + L*S). It then 
follows that (IHOj) will be automatically satisfied, while (IHTj) then only requires that (700 = 
— (|L*L + zilf) in order to obtain the Lindblad generator C (X). 

This leads us precisely to the coefficients Gi{t) stated in the lemma, and the identity 
Tt{A ® X) + Qt{A ® X)iy^ (t) (a+ 8) /) + u*^ (t)* (a_ ® /) Qi(A ® X) = 6^^ (A ® X) (82) 

The following Bayes' relation is proven along similar lines to pj Theorem 6.2]. 

Lemma 3.10 For a^ 7^ 0, define 

Q(yl ® X) = (/ ® U*{t))Ej,[F*{t){A ® X)F(t) | / ® ^t](/ ® t/(t)). (83) 
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In order to determine the differential equation for ^^[A^X), we first define, for A0X & 
B(C^ (g) ijs) tlie process 

jt (A ® X) = E^[F{tyA ^XF{t)\I ® ^t], (85) 

so that q(A (g)X) = {I(g) [/(t))*7i(A (g)X){I(g) U{t)). We then have 

Lemma 3.11 Let a^ ^ 0. The process 7^ {A ® X) satisfies the QSDE 

d^t {A®X) = Tt [A ®X)dt + /3t {A ® X) dZ (t) (86) 

where Tt {A X) , f3t {A ^ X) G J^, are given by 

(3t{A^X) =7,(Q,(A®X)) 

+z/e(t)7t (^ ® ^ (^+ ® /)) + 1^*^(0*7* ((^- ® /) A ® X) - 7t (A ® X) 0i, 
n(A®X) = 7i(Ti(A®X)) 

+<{tht {Qt {A ® X) a+) + z/*e(t)*7* (^-Qt (^ ® ^)) - A (A ® ^) ^t, 

with 

Ot = E^[ijyCit)a^ + Jy*Cit)cr-) ® I\I ® ^t]. 

Proof. Setting Rt = F{t)*{A ® X)F(t), we have that 

dRt = F{t)*Tt {A (g) X) F{t)dt + F{t)*Qt{A O X)F{t)dZ (t) 
and our aim is compute 74 (A (g) X) = Ex;[i?i|/ ® ^]. In particular, 

E^[{Rt--ftiA(S)X))Dt]=Q (87) 

for every Dt E I ® !^t and we now apply a technique similar to the characteristic function 
method, this time using the input process Z and taking the process Dt to satisfy the QSDE 
dDt = f (t) DfdZ (t) with Dq = I, for given integrable /. From the Ito product rule we 
then have 

= Es [{dRt - d^t {A ® X)) Dt + (Rt - 7* (^ ® X)) dDt + {dRt - d^t (A ® X)) dDt] 

(88) 
and making the ansatz that d'jt {A ® X) = Tt{A^ X) dt + l3t{A0 X) dZ (t) for unknown 
coefficients Tt {A X) and (3t {A ® X) we see that 

= Es[(F(t)*Ti {A ® X) F{t) -Tt{A^ X))Dtdt 

+ {F{tyQt {A ® X) F{t) - A (A ® X))DtdZ (t)] 
+E^ [{Rt - 7t (^ ® ^)) A/ (t) dZ{t)] 
+E^ [[F(t)*Qi (A ® X) F{t) -^t{A® X)] A/ {t) dt {t)] . 

We now make use of Lemma 13.31 again and apply the commutation relations F{t)a^ = 
a+F{t), a^F{t)* = F*{t)a^. (Note that a+ will not commute with F* {t).) Inserting 
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Es [-1/ ® ^t] under the expectation sign, then separating coefficients of Dt and Dtf{t), we 
obtain the equations 

= 7t (Tt (A ® X)) -Tt{A^ X) 

+uat {Qt {A ® X) {a+ ® /)) + u^t ((a_ ® /) Qt (A^X)) 
-(5t {A ® X) Es H{t)a^ ®I + u*i{tya^ ® I\I ® S't] , 

= ^t{Qt{A®X))-l5t{A®X) 

+vi{t)-it {A®X (a+ ® I)) + v*i{ty-it ((a_ ®I)A®X) 
-7i (A ® X) Es K(0^+ ® / + '^*eW V- ®I\I®^t]- 

Rearranging these expressions yields the relations in the statement of the lemma. 

Theorem 3.12 Let Uq ^ 0. The unnormalized conditional expectation q(A ® X) defined 
by ( fg^) satisfies the equation 



dqt {A®X)= Q {gt{A ® X)) dt + \t{A® X) dY{t), (89) 



where 



Xt {A®X) = q(A ® XLt + L;A ®X)-qt{A® X) Ku (90) 

Lt = I®L + ut^{t)a+ ® S, (91) 

dY{t) = dY (t) - Ktdt, y (0) = 0, (92) 

Kt = c^t {Hit)<y+ + ^*e(i) V_) ® I) . (93) 

Proof. We remark that by inspection the coefficients in the QSDE for 7^ {A ® X) 
simplify to 

Pt{A®X) = ^t{Qt{A®X)) 

+pmit {A®X (a+ ® I)) + z/*e(i)*7t ((^- ®I)A®X)--ftiA®X) Ot 
= jt{A ®X{I®L + Uta+®S) + {I®L* + i/*e(t) V_ ®S*)A® X) 
-7,(A®X)^„ 
Tt iA®X) = -ft (Tt (A ® X)) + u^{t)^t (Qt (A ® X) a+) 
+u*^{ty-ft {a^Qt {A ® X)) -f3t{A® X) Ot 
= ^t{g{A®X))-l3t{A®X)et, 

and therefore 

d^t {A®X)= 7i {g{A ® X)) dt + /3t{A® X) [dZ (t) - Otdt] . 

The QSDE for q(A ® X) = (I ® U{t)y-ft{A ® X){I ® t/(t)) is then readily deduced 
from the unitary rotation noting that Kt = (I ® U{ty*6t{I ® f/(t)) and Xt{A®X) = 

{i®u{ty*i3t{A®x){i®u{ty. 
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Corollary 3.13 For ao 7^ 0, the conditional expectation ntlA (g) X) defined by ([5^ and 
given by (84^ satisfies equation ^5^) derived in Theorem \3.6[ 

Proof. We see that c^q (/ (g) J) = A* (/ ® I) dY (t) and so 



d- 



1 



dY (t) + 



Q (/»/)■ 



rXt (J ® J)' dt. 



A, (/ 8) /) 
However, we note from (I5U1) that 

By an application of the Ito product rule, the normalized filter therefore satisfies 

dTTt{A (g) X) 
= 7rt{gt{A®X))dt 

I Q (/ 8) /) J L Q (/ ® /) _ 

-TTt (A X) TTt (Zt + L;) } [dY{t) - nt (it + L*) dt 
= ^i(^t(A X))dt + ?{i(^ ® X)rfM/ (t) , 

since we have from ( 19T|) . (l56l) and (!57|) 

?/i(A ® X) = TTi (a ® XLi + L*A ® X) - TTi (A ® X) Tit (it + ^i*) , 

dW{t) = dY -Tit (it + It) dt. 



Therefore we recover equation ( l55l) . 



4 Fields in a Superposition of Coherent States 

4.1 Superposition of Coherent States 

In this section we take the field to be in a superposition state 



(94) 



where \fj) are coherent states and the complex numbers Oj {j 
normalized weights (described further below). 



1 , . . . , n) are non-zero 
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Coherent vectors |/) may be expressed in terms of the vacuum vector using the Weyl 
(or displacement) operator [28] W{f) which serves as a "density": 

\f) = W{f)\0). (95) 

While the collection of all coherent vectors is dense in the Fock space, they are not orthog- 
onal, and indeed the inner product (in the Fock space) is given by 

{f\9) = eM-l\\f\\l-l\\9\\l+{f,9)2). (96) 

Here, || / \\l= (/, /)2 and (/, g)2 are the -^^([0, oo), C) norm and inner product respectively. 
The superposition state \ip) given by flMl) is specified by a choice of coherent vectors \fj), 
with weights aj ensuring normalization: (ipli/j) = Ylijk^*j^kgjk = 1, where Qjk = {fj\fk)- 

For a system operator X acting on S)s, and F is a field operator acting on the Fock 
space ^, the expectation with respect to the state |?7) ® I^E') is defined by 

Er,^[X®F] = {r]^\{X ^ F)\r]^) = {r]\X\r]){^\F\^) 
= {7^\X\v)J2^*Mfj\F\fk) 

jk 

= J2^*j^k^jk[X^F], (97) 

jk 

where 

E^k[X®F] = {r]\X\rj){f^\F\f,) (98) 

for J, A; = 1, . . . , n. We write ]Eoo[-^ ^ F] = {r]\X\ri){0\F\0) for the vacuum case. 

Consider now the expectation of an adapted operator K{t) on the composite system 
S) = {S^s ® !?<]) ® d{t', this means that K{t) acts trivially on the future component ^(^t- 
Let X[o,t] is the indicator function for the time interval [0,t]. Now coherent vectors and 
Weyl operators factorize as |/) = \fX[o,t]) ® |/X(t,oo)) and W{f) = W{fx[o,t])'S)W{fx{t,oo)), 
respectively. Write 

W,-if) = Wifxm), W,+ {f) = Wifx(t,oo)). (99) 

Then we can express the coherent expectations of adapted processes K{t) in terms of the 
vacuum: 

E,,[Kit)] = Eoo[Wrifj)Kit)Wiifk)y\t) (100) 

where r^''{t) = (0|Vr+*(/j)Vr+(/fc)|0) satisfies 

r^'it) = -{f;{t)fk{t) - llfM' - l\fk{t)\y{t), r^\0) = 1. (101) 

Note that j = k is the standard coherent expectation, in which case r^^{t) = 1. 

The following lemma shows how expectations of stochastic integrals with respect to the 
superposition state can be evaluated. 
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Lemma 4.1 Let K{t) he a hounded quantum stochastic process defined hy [T^) . where Mq, 
M± and Mi are hounded and adapted. Then we have 

Ejk[K{t)] = E,k[[ Mo{s)ds]+ [ M4s)Ms)ds 
Jo Jo 

+ [ M^{s)f;{s)ds+ [ Mi{s)f;{s)fk{s)ds]. (102) 



Proof. Equation (11021) follows from the following eigenstate property of coherent 
vectors: 

dB{t)\f) = f{t)\f)dt, 

dA{t)\f) = dB*{t)f{t)\f). (103) 

D 



4.2 Embedding 

For the superposition of n coherent states, we use an n-level ancilla system, leading to the 
extended space 

^ = C"(g)io = i^©io©---©io (n times). (104) 

As in the single photon case, we allow the extended system to evolve unitarily according 
to / ® U{t), where U(t) is the unitary operator for the system and field, given by the 
Schrodinger equation ([2]). Let \ej), j = 1, . . . , n, be an orthonormal basis for C". We 
initialize the extended system in the state 

\a\ ^ — ' 
j 

where aj ^ for all j and |ap = ^ • a*aj (so that (S|S) = 1). This state evolves according 
to |S(t)) = (/©[/(t))|S). 

Let A be an operator acting on C", i.e. a complex nx n matrix, A = (ajk), j,k = 
1, . . . ,n. Then expectation in the extended system is defined by 

Es[^ © X © F] = (S|(A © X ® F)|S) = —J2ajka*akEjk[X © F]. (106) 

jk 

Expectations of quantum stochastic integrals can be compactly expressed in the ex- 
tended system, as the following lemma shows. 
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Lemma 4.2 Let M{t) be adapted. Then 

Es[ / A® M{s)dB{s)] = Es[ / (AC(s)) ® M{s)ds], (107) 

Jo Jo 

Es[ / ^ ® M(s)dS*(s)] = Es[ / (C^"(s)A) ® M(s)(is], (108) 

Jo Jo 

Ee[ / A ® M(s)rfA(s)] = Es[ / iC\s)AC{s)) M{s)ds], (109) 

Jo Jo 

where 

C(t) = diag[/i(t),...,/„(t)]. (110) 

Notice that the expectations of the stochastic integrals are expressed in terms of the 
action of the matrix C{t) on the ancilla factor A. 

4.3 Master Equation 

In this section we show how the the unconditional expectation 

^it{X) = E,^[X{t)] (111) 

may be computed from a collection of differential equations. We do this through a differ- 
ential equation for the unconditional expectation 

lit{A®X)=E^[A(^X{t)] (112) 

for the extended system. 

Let R he an n X n matrix defined by Rjk = 1 for all j,k = 1, . . . ,n, and define 

GtiA 0X) = A® C{X) + {AC{t)) ® [L\ X]S + {C\t)A) ® S*[X, L] 

+{C\t)AC{t))®{S*XS-X)). (113) 

Lemma 4.3 The unconditional expectation lilll\) with respect to the superposition state 
|\E') (defined by [94\)) is given by 

and the master equation for the expectation I^11S\} in the extended system is 

^MA ® X) = fitiGtiA X)), (115) 

at 

with initial condition fLo{A ® X) = j^{r]\X\ri) J^jk^jk'^'j'^k- 
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Proof. By definitions fll06p and (|97]) we have 

Es[i?®X(t)] = ^^a*afcE,4X(t)] (116) 

jk 

^ E,*[X(t)], (117) 



laP 



and in particular 

Es[i?®/] = ^. (118) 

From these expressions, we see that 

E,^[X{t)] = \a\^E^[R ® X{t)] = ^i^^^^f^^ , (119) 

which proves (I114p . 

The differential equation (IllSp follows from the QSDE ([3]) for X{t) = jt{X) and rela- 
tions (ITU71) -(ITUgD upon evaluating the differential dEj^[A ® X{t)]. D 

Theorem 4.4 The unconditional expectation iJ,t{X) when the field is in the superposition 
state |\E') (defined by (94\)) is given by 






MX) = -^^ \ '\:,: . (120) 



where fil {X) is given by the system of equations 

j^f4\x) = f4\Qt{x)), (121) 

and where 

GtiX) = C{X) + S*[X, L]f;it) + [L\ X]Sfkit) + iS*XS - X)f*{t)Ut). (122) 

The initial conditions are 

4\X) = {v\X\rj)gj,. (123) 

Proof. Define 

fitiX) = E,,[Xit)]. (124) 

Then as in the proof of Lemma 14.31 we may show that 

fii\X) = ^f,,i\e,){e,\^X). (125) 

The the relation 012Up follows from 0114p . The differential equation 01211) follows from 
equation (IllSp with A = \ej){ek\- □ 
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4.4 Superposition State Filter 

In this section we show how the conditional expectation X(t) = T^ti^) defined by ([7]) can 
be evaluated using a system of conditional equations. This will make use of the conditional 
expectation 

ntiA®X)=Ej:[A^X{t)\I^^t]. (126) 

for the extended system. 

Lemma 4.5 The conditional expectation X{t) = TitiX) defined by ^ with respect to the 
superposition state |\1/) is given by 

,,(X) ^ M|^. (127) 

The quantum filter for the conditional expectation 7tt{A X) is 

dnt{A 0X) = ntigt{A ® X))dt + Ht{A ® X)dW{t) (128) 



with initial condition t^q{A ® X) = j^{ri,Xr]) Yljk'^jk(^*j^k, wh 



ere 



nt{A ®X) = ^t{A ®X{I(^L + C{t) ®S) + {I®L* + C\t) ® S*)A ® X) 

-nt{A ® X)nt{I ®L + C{t) ® S + I(^ L* + C^{t) (g) S*) (129) 

and W{t) is a '3^t- Wiener process given by 

dW{t) = dY{t) - nt{I ®L + C{t) ®S + I®L* + C^{t) O S*)dt, W{0) = 0. (130) 

The initial condition is tto{A (g) X) = j^{ri\X\r]) J^jk'^jk'^^^k- 

Proof. Let K e 'S^t- Then we have 

E^[nt{R^X){I^K)] = Es[(i?®X(t))(/®K)] 

= T^y,(^;ak^jk[X{t)K] 

jk 
^ Er^^lX {t)K] 



lap 
-^E,^[nt{X)K] 



\a\ 
Ej:[R(^Trt{X)K] 

Es[Ej:[R ® irtiX)K\I ^ -^t\ 
Es[^i(i?®/)(/®7rt(X))(/®ir)]. (131) 



This proves (11271) . 
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The filtering equation f ll28p is derived using the characteristic function method [29] , 
[7], [5]. We postulate that the filter has the form 

dnt{A ® X) = Tt{A ® X)dt + Ht{A ® X)I ® dY{t), (132) 

where J-^ and Tit are to be determined. 

Let / be square integrable, and define a process c/(t) by 

dcfit) = fit)cf{t)dY{t), Cf{0) = 1. (133) 

Then I ®Cf{t) is adapted to I®'3^ti and the definition of quantum conditional expectation 
[9|, sec. 3.3] implies that 

Es[A ® (X(t)c/(t))] = Es[fri(A ® X{t)){I ® c/(t))] (134) 

holds for all /. By calculating the differentials of both sides, taking expectations and 
conditioning we obtain 

Es[A ® {dX{t)cf{t))] = Es[(/ ® Cf{t))%t{G{A ® X)) (135) 

+ (/ ® f{t)cf{t))nti{A ® X)(/ ® L + Cit) ® 5) 
+ (/ ® L + Ct(t) ® 5*)(A ® X))]dt 

and 

Es[A®(d^i(A®X)c/(t))] (136) 

= Es[(/ ® c/(t){j;(A ® X) + Hi(A ® X)^i(/ ® L + C(t) ® S + / ® L* + C\t) ® 5*)} 

+ (/ ® f{t)cf{t)){ nt{A ® X)^i(/ ® L + C{t) ®S + I^L* + C\t) ® 5*) + UtiA ® X)}]rfi. 

Now equating coefficients of Cf(t) and f{t)cf{t) we solve for J-i(A ® X) and 'Hi(A X) to 
obtain the filtering equation 



We now show that W{t) is a ^-martingale, and since dW{t)dW{t) = dt, then by Levy's 
theorem [12] we have that W{t) is a ^-Wiener process. Indeed, for any K G ^ we have 

E4{I®dW{t)){I®K)] 
= Es[(/ ® rfr(t) - ^t(/ ® L + C{t) ^S + I^L* + C\t) ® 5*)rft)(/ ® ir)] 
= Es[/ ® (L(t) + L*(t)) + C(t) ® 5 + C*(t) ® 5* 

-^i(/ (g) L + C(t) (g) 5 + / ® L* + C^(t) ® 5*)dt)(/ (g) K)]dt = 0. (137) 

This completes the proof. D 

Theorem 4.6 The unconditional expectation fit{X) when the field is in the superposition 
state |\1') (defined by [94\)) is given by 



E,.c.*«,7rf(X) 
E,,«*a,7rf(/) 



MX) ^ t:' r z:: ^ (m 
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where the conditional quantities vr^ (X) are given by 

rfvrf (X) = 7ri\g^\X))dt + (vrf (X(L + SMt)) + (L* + S*f;{t))X) 



J 



r/ie innovations process W{t) is a '2Vt Wiener process with respect to the superposition state 
1^) and is given by 



laJ^ 



dW{t) = dY{t) - J2 TjV^t'(^ + '^/i(^) + ^* + S*fjii))dt- (139) 

j 

The initial conditions are 

<\X) = {v\X\rj)g,,. (140) 

Proof. These assertions follow upon substitution of 

ntiX) = ^n,{e,el^X). (141) 

into the relevant expressions from Lemma I4.5[ D 

5 Discussion and Conclusion 

In this paper we have derived the master equation and quantum filter for a class of open 
quantum systems that are coupled to continuous-mode fields in non-classical states: (i) 
single photon states, and (ii) superpositions of coherent states. The quantum filter in both 
of the cases we consider consists of coupled equations that determine the evolution of the 
conditional state of the system under continuous (weak) measurement performed on the 
output field, in contrast to the familiar single filtering equation for open Markov quantum 
systems that are coupled to coherent boson fields. This coupled equations structure of the 
master and filter equations is a reflection of the non-Markov nature of systems coupled 
to the non-classical fields. Indeed, a key feature of our approach is the embedding of 
the system into a larger extended system, a technique often employed in the analysis 
of non-Markov systems, providing an elegant framework within which to study the the 
dynamics, both unconditional and conditional, of the system. In contrast to Markovian 
embeddings [10], [19], [20], the extended system (including the field) is initialized in a 
superposition state. This embedding provides a framework within which the tools of the 
quantum stochastic calculus may be efficiently applied to determine quantum filtering 
equations. We expect that the use of suitable embeddings, both Markovian and non- 
Markovian, could be adapted to study quantum systems that are coupled to other types 
of highly non-classical fields. 
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